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ABSTRACT 

We present the discovery of four remote star clusters in M33, one of which is 
of an extended nature. Three of the clusters were discovered using survey data 
from the Isaac Newton Telescope Wide-Field Camera while one was discovered 
serendipitously in a deep image taken with the Hubble Space Telescope's Ad- 
vanced Camera for Surveys. With projected radii of 38-113 arcmin (9.6-28.5 kpc 
for an assumed M33 distance of 870 kpc), these clusters lie significantly beyond 
all but one of the currently-confirmed clusters in M33. The clusters have mag- 
nitudes and colors consistent with their being old to intermediate-age globular 
clusters. Indeed, they bear a strong resemblance to the outer halo GC population 
of the Milky Way and M31 in terms (V-I)o colours. The three outermost clusters 
are projected on the far side of M33 with respect to M31, an asymmetry that 
could suggest tidal interactions have affected M33's globular cluster distribution 
at large radii. 

Subject headings: galaxies: globular clusters: - galaxies: individual (M33) - galaxies: 
formation - galaxies: evolution - galaxies: star clusters - galaxies: stellar content 
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Introduction 



Globular clusters (GCs) have proved to be valuable probes of galaxy formation, 
providing key clues about the assen ibly histories of the Milky Way (hereafter MW) (e.g. 



Searle fc Zinn 



19781), and M31 (e.g. 



Perrett et al. 



2OO2I ) . A particular property of GCs that 



makes them especially important is that they can often be found out to extremely large 
radii in galaxy halos. Such outer halo GCs are key tracers of the bar yonic content and 



dyna mics of galaxies on scales of tens to hundreds of kiloparsecs (e.g. 
2000). 



Evans fc Wilkinson 



First discovered by .Madore fc ArpI (119791 ). AM-1 remains the most distant GC known 



in the MW with a galactocentric distance of 123 kpc. M31 hosts GCs at much larger radial 



distances, out to > 200 kpc, however this population has only been u ncovered in t 



le last 



few years thanks to sensitive panoramic digital imaging surveys (e.g. 



Martin et al. 



2006 



Huxor et al. 



20081 ). For example, using 84 square degrees of imaging data from the Isaac 



Newton Telescope Wide-Field Came ra (INT/WFC) an d Canada France Hawaii Telescope 



(CFHT) MegaCam surveys of M31, 



Huxor et al. 



(I2OO8I ) were able to increase the number 



of known GCs with projected radii > 1 degree (~ 14 kpc at M31's distance) by more than 
75%; this is especially remarkable given that most of this data came from a single outer 
quadrant of M31 and suggests that many more halo GCs remain to be discovered in that 
system. 

M33 is the least massive spiral in the Local Group and there has been much interest in 
understanding how its star cluster system compares to that of its larger brethren. Previous 
studies have suggested that, like the MW and M31, M33 has a population of star clusters 
which exhibit halo kinematics but that these clusters possess a much larger age spread 



and a higher mean 



Chandar et al 



metallicity than their MW counterparts (e.g. 



Sarajedini et al. 



2OO0I : 



20021 ). Indeed, M33 appears to have a cluster population with a very broad 
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range of ages, bearing a stronger resemblance to the cluster populations of the lower mass 
Magellanic Clouds rather than the MW or M31. 

One drawback of all M33 star cluster studies to date is that they have focused only on 
objects projected on the inner regions of the galaxy, with little exploration of the galaxy 
outskirts. This bias has made it rather d i fficult to unambiguously associate clusters with 

(I2OO2I ) study, for example, includes only objects 



the disk or the halo. The 



Chandar et al. 



that lie within a projected radius of 5 kpc, c orresponding to just ove r half the isophotal R25 



radius of 9 kpc fide Vaucouleurs et al.lll991l) . iSarajedini fc Mancond (120071 . hereafter SM07) 



have recently compiled a list of all M33 star clusters and cluster candidates identified in all 
previous studies. This has led to the creation of a unified catalog of ~ 300 'high confidence' 

SM07 catalog lies at a mere 



Zloczewski et al. 



( 12008 



hereafter 



objects. The most remote cluster reported in the origina., 
7 kpc in projected radius. Since the publication of SM07, 
ZKH08) have conducted a search for new star clusters in the outskirts of M33 based on 
0.75 square degrees of deep CFHT Megacam imagery. Although their survey has led to the 
discovery of more than 120 new star cluster candidates, none of these lie beyond a projected 
galactocentric radius of 8.3 kpc. 

In order to establish whether M33, like its larger neighbours, also possesses a population 
of remote halo star clusters, we have undertaken a search for GCs in INT/WFC survey data 
of M33, as well as in our own deep Hubble Space Telescope Advanced Camera for Surveys 
(HST/ACS) imagery of two outer fields in M33. Our search has resulted in the identification 
of four new outlying star clusters in M33, which lie at projected radii of 10 -30 kpc. A fifth 



cluster was also uncovered in our search but it has since been published by 



Stonkute et al. 



(I2OO8I ) who independently found the object in Subaru Suprime-Cam data (M33-EC1). 
Throughout this paper, we adopt a distance modulus to M33 of (m — M)o = 24.69, 
corresponding to 870 kpc, for consistency with SM07. 
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2. Observations 

The primary dataset used in our GC search is the INT/WFC survey of M33 which was 
conducted over several observing runs in 2002-2008. The WFC is a four-chip CCD mosaic 
camera with a 0.29 deg^ field of view and a pixel scale of 0.33 arcsec pixel~\ equipped to 
the INT 2.5m telescope on La Palma. Exposures of 900-1200 seconds were taken in both 
the Johnson V and Gunn i bands, with a median seeing of 1.2". This depth is sufficient to 
detect individual stars in M33 to V~ 24.5 and i~ 23.5 with a typical signal-to-noise of 5. 
The fields observed comprise a contiguous, almost circular, region, covering approximately 
12 square degrees, very much greater in extent than any previously-published survey of 
M33 (see Figure [1]). The survey data were processed using the standard INT Wide Field 
Survey pipelin e prov ided by the Cambridge Astronomical Survey Unit; as discussed in 



Irwin fc Lewis! ( l200ll ). this pipeline provides basic data processing and astrometric and 



photometric calibration, as well as object detection and classification. 

Our star cluster search concentrated on only those regions beyond the main optical 
disk of M31. We adopted a strategy similar to that used in our search for GCs around 

I ir^ 

M31 fiHuxor et al.ll2008l ) and which exploits the fact that clusters at the distance of M31 
and M33 should be just resolved in good seeing. Indeed, given the survey depth, we 
expect detection of individual red giants down to ~ 3 magnitudes below the tip of the 
red giant branch. Initial GC candidate selection was based on source magnitude, color 
and morphological classification (stellar/non-stellar) provided by the photometric pipeline, 
exploiting the fact that most known GCs lie within specific ranges of absolute magnitude 
(-10.5 < Mv < -3.5, equivalent to 14 < F < 21 at the distance of M33) and color 
(0<y — /<1.7). These generous regions of parameter space are sufficient to allow for 
photometric errors and avoid exclusion of any slightly atypical GCs from the candidate list 
generated. Unlike our search for M31 GCs, ellipticity and object FWHM were not employed 
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as search criteria since the number of initial candidates selected using just magnitude, 
color and morph ology alone was not prohibitively large and all candidates could be visually 



inspected. As in 



Huxor et al. 



(l2008l ). every single survey image was also visually inspected 



in order to identify well-re solved diffuse extended GCs. Such objects are rather common in 



the outskirts of M31 (e.g. 



Huxor et al 



2005 



20081 ) and would not have been automatically 



found by the method above as they do not appear as a single source in the object catalog. 

We also conducted a search for outlying star clusters in two deep HST/ACS pointings 
that we obtained along M33's northern major axis at radii of 36' and 46'. These data were 
taken in the F606W and F814W filters as part of GO PID 9837 (PI = A. Ferguson) with 
total exposure times of 18050 and 12955 sees respectively and are being used for an analysis 
of the star formation history of M33's outer disk (Barker et al., in preparation). The star 
cluster search in this dataset was conducted purely by visual inspection of the drizzled 
stacked frames since star clusters at the distance of M33 are very well-resolved in HST/ACS 
imagery. 



3. Results 

Our search resulted in the identification of four new outlying star clusters in M33, 
three of which were discovered in the INT/WFC survey data and one in the HST/ACS 
data. It is worth noting that although the HST/ACS frames lie within the INT/WFC areal 
coverage, the HST cluster was not independently found in the INT/WFC survey data. This 
cluster is faint (see below) but, more importantly, rather compact in appearance hence it 
was not possible to distinguish it as a cluster in our ground-based imaging. Figure [H shows 
the locations of these objects with respect to the 'high confidence' clusters from SM07 and 
the ZHK candidate clusters. Our new discoveries lie at galactocentric radii ranging from 
38-113 arcmin, which corresponds to 9.6-28.5 kpc for an assumed M33 distance of 870 kpc. 
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and are clearly more remote than all but one of the known population. If these clusters are 
associated with M33's disk rather than halo, then their deprojected radii would be even 
more extreme, corresponding to 10-38 kpc. 

Figure [2] shows images of the four newly discovered clusters. Two of these clusters 
(HM33-B and HM33-C) are high surface brightness and compact while one clust er 



Stonkute et al 



fl2008h . 



(HM33-A) is diffuse and extended, much like M33-EC1 reported by 
Figure 2 also shows the stacked F606W image of the cluster (HM33-D) discovered in the 
HST/ACS image. 

Integrated photometry of the new INT/WFC clusters was undertaken with the IRAF 
task apphot 3. Fluxes within an aperture of radius 8" were obtained for clusters HM33-B 
and HM33- C; this apertu r e was selected to allow direct comparison with our photometry 



of M31 GCs 



Huxor et al. 



(120081 ). The (V-I) color for these clusters was measured within a 
smaller aperture of 4" radius that as chosen to reduce the error from the background and 
assumes that there is no color gradient in the clusters. For the extended cluster IIM33-A, 
the magnitudes were determined within a larger aperture of 12" radius so as to enclose 
the bulk of the light. The INT/WFC survey uses Johnson V (V) and Gunn i passbands 



hence a color trans 
V and I. Following 



brmation was requi r ed to obtain the Johnson/ Cousins equivalents 



McConnachie et al. 



(120031 ). we adopted V = V + 0.005(\/ - J) and 



I = i — 0.101(y — I). Due to the high luminosity of the clusters, the formal errors reported 
by apphot are quite small. The cluster magnitude uncertainties are dominated by the 
zero-point errors of ~0.02 mag while colour errors are estimated to be ~ 0.03 mag. These 
errors do not account for contamination from foreground stars or background galaxies as 



^IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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these are difficult to estimate. Finally, the cluster magnitud es and colors were corrected for 



extinction by interpolating within the 



Schlegel et al. 



(119981 ) mapqj to the position of each 



cluster. The INT/WFC clusters have absolute magnitudes in the range M^o ~ —6 to —7 
placing them near the peak of the GC luminosity function. 



For the HST/ACS cluster, we followed a similar procedure to that of 



Sarajedini et al. 



(120061 ). First, the F606W drizzled image was convolved with a Gaussian kernel (1 a 



20 pixels) as otherwise the light distribution was too patchy to determine the cluster 
centre by eye or by fitting a Gaussian. Once the center was determined, magnitudes were 
measured on two of the drizzled sub-images produced by the STScI CALACS pipeline, with 
total exposures of 5240 and 10480 seconds for F606W and F814W filters respectively. We 
determine fiuxes in a circular aperture of radius 3.5"; as the cluster lies near the edge of an 
ACS chip, this was the ma ximum aperture size that could be chosen which lay entirely on 



the chip. Interpolating the 



Schlegel et al 



(119981 ) dust maps to the cluster's position gave a 



reddening of E(B-V) = 0.046. If we adopt the Sirianni et al. (2005) extinction ratios for 
a G2 star in the ACS/WFC filter system and the most up-to-date ACS zeropointjf] then 
the dereddened magnitude and color are F6O6W0 = 20.46 and (F606W-F814W)o = 0.77. 
The S/N of these images is very high and the formal uncertainties on the magnitudes are 
very small. To convert the ACS/WFC magnitudes to the Johnson- Cousins system, we 
used the synthetic transformations of Sirianni et al. (2005) which gave Vq = 20.60 (with 
a systematic zero-point error of 0.05 mag) and (V-I)o = 0.92 mags. This corresponds to 
a cluster magnitude of My^ of -4.09, the faintest yet known star cluster in M33, and is 
comparable to the low-luminosity MW cluster Pal 13. Table [H lists the positional data and 
corrected photometry for all clusters. 



^http: / / astro.berkeley.edu/ ~marc/dust / data/data.html 
^http:/ /www. stsci.edu/hst/acs/analysis/zeropoints. 
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4. Discussion 

In the following, the newly-discovered clusters are compared with the known cluster 
population in M33. For this purpose, we use a subset of 167 of the 296 'high confidence' 

magnitudes availab le. We corrected 



Schlegel et al. 



(Il998l ) dust maps to 



clusters in the SM07 catalogue that have both V and 
this sample for extinction by interpolating within the 
the position of each cluster. This resulted in a range of Ay from 0.14 to 0.66. We do not 
use the ZKH08 sample in these comparisons as the authors note that the new GCs in their 
catalog are only "probable", and may contain contaminants. 

Figure [3] shows the location of the new star clusters on a plot of Vq magnitude versus 
(V— l)o color. There is a broad spread in magnitudes at all colors and the new objects 
sit comfortably within this distribution. On the other hand, the HST/ACS cluster is far 
fainter than any previously-confirmed M33 star cluster. However, the color of HM33-D, 
combined with its compact appearance on the HST/ACS images, supports its identification 
as an intermediate or old age GC. This conclusion is further strengthened by the isochrone 
fits to the color-magnitude diagram which indicate an age of > 3 Gyr and [M/H]> —0.7 
dex (Barker et al., in preparation). The existence of such a cluster supports the notion that 
a population of very faint GCs may await discovery in deeper panoramic imaging of the 
outskirts of M33. 



Stonkute et al. 



fl2008h . show a 



Our newly-discovered clusters, along with that of 
notable lack of spread in (V— l)o color when compared to the remaining M33 cluster 
population. This is particularly apparent when (V— l)o color is plotted against projected 
radial distance from M33 (Figure H]). The new clusters are slightly redder than the overall 
cluster population having an average (V— I)o color of 0.88 ± 0.05 compared to the average 
(V-I)o for the inner population being 0.67±0.30 mags. The large spread in color at smaller 
radii is not surprising since the SM07 sample contains many young clusters, which have 
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bluer colors than typical old GCs, and reddening issues will be more problematic here. The 
colors of the five outermost clusters in Figure 4 are very similar to those of M31 halo GCs 
with projected radii greater than 30 kpc, which have a mean value for (V-I)o = 0.91 ± 0.15, 
as well to Milky Way clusters beyond 15 kpc, which have (V-I)o = 0.92 ± 0.12 (Huxor et al. 
2009). That the MW and M31 outer halo GCs are gener ally old (> 10 Gyr) and metal-poor 



([Fe/H]< -1.5) fflarris 



1996 



Mackey et al. 



2006 



20071 ) supports the idea that the new 



M33 halo clusters are archetypical GCs t oo, albeit with c l uster HM33-A and M33-EC1 



having more extended structures. Indeed, 



Stonkute et al. 



(120081 ) were able to place some 



constraints on the age and metallicity of M33-EC1 from analysis of their color-magnitude 
diagram, finding an age > 7 Gyr and [M/H]< —1.4. 

Finally, it is notable that so few luminous objects have been discovered in the ~ 12 
square degree area we have surveyed around M33 (Figured]). Ignoring the central square 
degree of our survey which was not searched for new objects due to the presence of the 
bright disk, we calculate a GC surface density of ~ 0.4 deg~^ in the outskirts of M33 
which is almost half that derived for the GC surface density in M31 over the radial range 
30 — 100 kpc (Huxor et al. 2009). Given the extent of the M33 survey area, we can be 
confident that we have found most of the luminous GCs in the outer halo of the galaxy, 
although it is possible that a population of small, compact and faint clusters, such as the 
example found in our HST/ACS imagery, still remains to be uncovered. It is also rather 
interesting that the outermost clusters in M33 are all projected on the far side of the galaxy 
with respect to M31, despite the survey coverage being essentially uniform. Assuming the 
outermost clusters are moving on circular orbits with velocities of 80 km s~^, typical of 
the old halo population identified in the inner regions of M33 (Chandar et al. 2002), their 
orbital times will be ~ 1 — 2 x 10^ yr raising the question as to whether their present spatial 
configuration can be long-lived or just a chance coincidence. The clear asymmetry in the 
outer cluster distribution could be a signature that tidal stripping by its larger neighbour 
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may have dramatically affected the globular cluster distribution in M33. On the other 
hand, this may just be a chance alignment. Line of sight distances, radial velocities and 
sensitive cluster searches at larger radii will be required to test these ideas further. 

5. Summciry 

In this paper we have presented the discovery of four new outlying star clusters in 
M33 using INT/WFC survey data and deep HST/ACS imagery. Three of these clusters are 
compact while one is of the extended class. These objects, which have projected radii of 
38-113 arcmin or 9.6-28.5 kpc (for an assumed M33 distance of 870 kpc), lie significantly 
beyond all but one of the currently-confirmed star clusters in M33. Their colors and 
magnitudes are typical of the intermediate- age to old M33 GC population, as well as old 
age GCs in the far outskirts of the MW and M31. The surface density of halo GCs beyond 
the optical disk of M33 is roughly one half of that measured in M31 in the radial range 
~ 30 — 100 kpc. The new clusters have a strikingly asymmetric distribution around M33, all 
lying on the far side of the galaxy with respect to M31. Further observations are required 
to properly characterize the stellar populations of these new objects and determine their 
line-of-sight distances and kinematics. 

APH, AMNF and MKB arc supported by a Marie Curie Excellence Grant from the 
European Commission under contract MCEXT-CT-2005-025869. NRT acknowledges a 
STFC Senior Research Fellowship. We thank Dougal Mackey for useful discussions, and 
the anonymous referee whose comments greatly improved the paper. The Isaac Newton 
Telescope is operated on the island of La Palma by the Isaac Newton Group in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. This 
work is partly based on observations made with the NASA/ESA Hubble Space Telescope 
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Table 1. Basic Data for the New Outlying GCs 



ID 


RA 


Dec 


Vo 


(V-I)o 




Projected Galactocentric 


Deprojected Galactocentric 




J2000 


J2000 


mags 


mags 


mags 


Radius (kpc) 


Radius (kpc)* 


HM33-A 


01 35 41.78 


+28 49 15.5 


18.60 


0.83 


0.26 


28.6 


37.7 


HM33-B 


01 36 02.12 


+29 57 49.4 


17.59 


0.89 


0.17 


12.8 


20.4 


HM33-C 


01 37 14.46 


+31 04 27.8 


18.17 


0.93 


0.18 


12.7 


17.1 


HM33-D 


01 35 02.20 


+31 14 21.3 


20.60 


0.92 


0.13'' 


9.6 


9.6 



^Deprojected values calculated assuming j = 56 degrees and PA = 23 degrees. 
''Extinction in F606W. 



-1 - 



-2 - 




M33-EC1 



10 kpc ^ 



*A r - 1 



-1 



-2 



^ (degrees) 



Fig. 1. — The locations of the new clusters (black stars) are shown in relation to the 'high 



confidence' clusters in the SM07 cata. 
(blue circles). The 



Stonkute et al. 



ogue (red circles) and the 'probable clusters' of ZKH08 



(120081 ) cluster, M33-EC1, is labeled, and an arrow in- 



dicates the direction to M31. The solid ellipse indicates the R25 isophote and the dashed 



outline indicates the extent of our INT/WFC sur vey coverage 



clusters, along with the remote extended cluster of 
side of M33 with respect to M31. 



Stonkute et al. 



le outermost of the new 



(l2008l ). project on the far 



-17- 




Fig. 2. — Thumbnails (i-band) of the outer halo clusters found in the INT-WFC survey (A, 
B and C) and in the F606W HST/ACS data (D). The INT images span 30" x 30" while the 
ACS image spans 7" x 7". North is up and east is to the left. 
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Fig. 3. — Plot of extinction-corrected color versus apparent magnitude for the new INT/WFC 
clusters (inverted triangles), the new ACS cluster (triangle) and the 'high- confidence' SM07 
sample (open diamonds). 
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Fig. 4. — Plot of the extinction-corrected color versus radius for the cluster population in 
M33. The dashed line indicates the mean extinction-corrected color of the outer halo globular 
clusters in the MW and M31 from Huxor et al. (2009). 



